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Abstract 

Cold dark matter may be made of superweakly-interacting massive particles, super WIMPs, that 
naturally inherit the desired relic density from late decays of metastable WIMPs. Well-motivated 
examples are weak-scale gravitinos in supergravity and Kaluza-Klein gravitons from extra dimen- 
sions. These particles are impossible to detect in all dark matter experiments. We hnd, however, 
that superWIMP dark matter may be discovered through cosmological signatures from the early 
universe. In particular, superWIMP dark matter has observable consequences for Big Bang nu- 
cleosynthesis and the cosmic microwave background (CMB), and may explain the observed under- 
abundance of 7 Li without upsetting the concordance between deuterium and CMB baryometers. 
We discuss implications for future probes of CMB black body distortions and collider searches for 
new particles. In the course of this study, we also present a model-independent analysis of entropy 
production from late-decaying particles in light of WMAP data. 

PACS numbers: 95.35.+d, 98.80.Cq, 26.35.+C, 98.80.Es 
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I. INTRODUCTION 



Recently, we proposed that dark matter is made of superweakly-interacting massive par- 
ticles (superWIMPs) Q. This possibility is realized in well-studied frameworks for new 
particle physics, such as those with weak-scale supersymmetry or extra spacetime dimen- 
sions, and provides a qualitatively new possibility for non-baryonic cold dark matter. 

The basic idea is as follows. Taking the supersymmetric case for concreteness, consider 
models with high-scale supersymmetry-breaking (supergravity models) and i?-parity con- 
servation. If the lightest supersymmetric particle (LSP) is the neutralino, with mass and 
interaction cross section set by the weak scale M wca k ~ 100 GeV — 1 TeV, such models are 
well-known to provide an excellent dark matter candidate, which naturally freezes out with 
the desired relic density (H |j| . 

This scenario relies on the (often implicit) assumption that the gravitino is heavier than 
the lightest standard model superpartner. However,even in simple and constrained super- 
gravity models, such as minimal supergravity [H, H, El, 0] , the gravitino mass is known only 
to be of the order of M weak and is otherwise unspecified. Given this uncertainty, assume 
that the LSP is not a standard model superpartner, but the gravitino. The lightest standard 
model superpartner is then the next-lightest supersymmetric particle (NLSP). If the universe 
is reheated to a temperature below ~ 10 10 GeV after inflation (if, the number of gravitinos 
is negligible after reheating. Then, because the gravitino couples only gravitationally with 
all interactions suppressed by the Planck scale Mp\ ~ 1.2 x 10 19 GeV, it plays no role in 
the thermodynamics of the early universe. The NLSP therefore freezes out as usual; if it 
is weakly-interacting, its relic density will again be near the desired value. However, much 
later, after 

M * 1 ~10 5 s-10 8 s, (1) 



^weak 



the WIMP decays to the LSP, converting much of its energy density to gravitinos. Gravitino 
LSPs therefore form a significant relic component of our universe, with a relic abundance 
naturally in the desired range near Odm — 0.23 [9J. Models with weak-scale extra dimensions 
also provide a similar dark matter particle in the form of Kaluza-Klein gravitons 0, with 
Kaluza-Klein gauge bosons or leptons playing the role of WIMP [lOj . As such dark matter 
candidates naturally preserve the WIMP relic abundance, but have interactions that are 
weaker than weak, we refer to the whole class of such particles as "superWIMPs." 

WIMP decays produce superWIMPs and also release energy in standard model parti- 
cles. It is important to check that such decays are not excluded by current constraints. The 
properties of these late decays are determined by what particle is the WIMP and two param- 
eters: the WIMP and super WIMP masses, mwiMP and mswiMP- Late-decaying particles in 
early universe cosmology have been considered in numerous studies [ll, 12[ 3, IJ, EH EH Ell- 
For a range of natural weak-scale values of m-wiMP and ^swimp, we found that WIMP — > 
SWIMP decays do not violate the most stringent existing constraints from Big Bang nu- 
cleosynthesis (BBN) and the Cosmic Microwave Background (CMB) [1]. SuperWIMP dark 
matter therefore provides a new and viable dark matter possibility in some of the leading 
candidate frameworks for new physics. 

SuperWIMP dark matter differs markedly from other known candidates with only grav- 
itational interactions. Previous examples include ~ keV gravitinos [18], which form warm 
dark matter. The masses of such gravitinos are determined by a new scale intermediate 
between the weak and Planck scales at which supersymmetry is broken. Superheavy can- 
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didates have also been proposed, where the dark matter candidate's mass is itself at some 
intermediate scale between the weak and Planck scales, as in the case of wimpzillas (l9| . 
In these and other scenarios |2(|, the dark matter abundance is dominantly generated by 
gravitational interactions at very large temperatures. In contrast to these, the properties 
of superWIMP dark matter are determined by only the known mass scales M wea k and M-p\- 
SuperWIMP dark matter is therefore found in minimal extensions of the standard model, 
and superWIMP scenarios are therefore highly predictive, and, as we shall see, testable. 
In addition, superWIMP dark matter inherits its relic density from WIMP thermal relic 
abundances, and so is in the desired range. SuperWIMP dark matter therefore preserves 
the main quantitative virtue of conventional WIMPs, naturally connecting the electroweak 
scale to the observed relic density. 

Here we explore the signals of superWIMP dark matter. Because superWIMPs have in- 
teractions suppressed by M P1 , one might expect that they are impossible to detect. In fact, 
they are impossible to detect in all conventional direct and indirect dark matter searches. 
However, we find signatures through probes of the early universe. Although the super- 
WIMP dark matter scenario passes present constraints, BBN and CMB observations do 
exclude some of the a priori interesting parameter space with mwiMP , "^swimp ~ ^weak- 
There may therefore be observable consequences for parameters near the boundary of the 
excluded region. Certainly, given expected future advances in the precision of BBN and 
CMB data, some superWIMP dark matter scenarios imply testable predictions for upcom- 
ing observations. 

Even more tantalizing, present data may already show evidence for this scenario. Late 
decays of WIMPs to superWIMPs occur between the times of BBN and decoupling. They 
may therefore alter the inferred values of baryon density from BBN and CMB measurements 
by (1) destroying and creating light elements or (2) creating entropy 2l]|. We find that the 



second effect is negligible, but the first may be significant. At present, the most serious 
disagreement between observed and predicted light element abundances is in 7 Li, which is 
underabundant in all precise observations to date. As we will show below, the superWIMP 
scenario naturally predicts WIMP decay times and electromagnetic energy releases within 
an order of magnitude of r « 3 x 10 6 s and Cem = ^em^wimp ~ 10~ 9 GeV, respectively. 
This unique combination of values results in the destruction of 7 Li without disrupting the 
remarkable agreement between deuterium and CMB baryon density determinations |l7| . 

We then discuss what additional implications the superWIMP scenario may have for 
cosmology and particle physics. For cosmology, we find that, if 7 Li is in fact being destroyed 
by WIMP decays, bounds on // distortions of the Planckian CMB spectrum are already 
near the required sensitivity, and future improvements may provide evidence for late decays 
to superWIMPs. For particle physics, the superWIMP explanation of dark matter favors 
certain WIMP and superWIMP masses, and we discuss these implications. 



II. SUPERWIMP PROPERTIES 



As outlined above, superWIMP dark matter is produced in decays WIMP — > SWIMP + S 1 , 
where S denotes one or more standard model particles. The superWIMP is essentially 
invisible, and so the observable consequences rely on finding signals of S production in the 
early universe. In principle, the strength of these signals depend on what S is and its initial 
energy distribution. For the parameters of greatest interest here, however, S quickly initiates 
electromagnetic or hadronic cascades. As a result, the observable consequences depend only 
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on the WIMP's lifetime r and the average total electromagnetic or hadronic energy released 
in WIMP decay [H Q EE Q IH, EE O H3 

We will determine r as a function of the two relevant free parameters mwiMP and mswiMP 
for various WIMP candidates. These calculations are, of course, in agreement with the 
estimate of Eq. JI}, and so WIMPs decay on time scales of the order of a year, when the 
universe is radiation-dominated and only neutrinos and photons are relativistic. In terms of 
r, WIMPs decay at redshift 

"10 6 



4.9 x 10 



T 



and temperature 



T 



90M? 



47r 2 T 2 g*(T) 



0.94 keV 



10 6 s 



where M* = M P \/\/87r ~ 2.4 x 10 18 GeV is the reduced Planck mass, and g*(T) 
the effective number of relativistic degrees of freedom during WIMP decay. 
The electromagnetic energy release is conveniently written in terms of 

Cem = ^em^wimp , 



(2) 

(3) 
29/4 is 

(4) 



where Eem is the initial electromagnetic energy released in each WIMP decay, and Iwimp = 
^wimp/^° is the number density of WIMPs before they decay, normalized to the number 
density of background photons n^ G = 2((3)T 3 /tt 2 . We define hadronic energy release simi- 
larly as Chad = ^had^wiMP- In the superWIMP scenario, WIMP velocities are negligible when 
they decay. We will be concerned mainly with the case where S is a single nearly massless 
particle, and so we define 

2 2 

TTi _ m WIMP — m SWIMP fr\ 

zmwiMP 

the potentially visible energy in such cases. We will determine what fraction of Es appears 
as electromagnetic energy e EM and hadronic energy e had in various scenarios below. For 
Yvvimp, each WIMP decay produces one superWIMP, and so the WIMP abundance may be 
expressed in terms of the present superWIMP abundance through 



^WIMP — ^SWIMP,t — 

~ 3.0 x 10~ 12 



^SWIMP, 

TeV 

^SWIMP 



^SWIMPPc 
^SWIMP^^ 

a 



SWIMP 



0.23 



(6) 



For Eem ~ Es ~ "^swimp ~ ^weak, Eqs. (J3J) and © imply that energy releases in the 
superWIMP dark matter scenario are naturally of the order of 



c 



EM 



10~ 9 GeV . 



(7) 



We now consider various possibilities, beginning with the supersymmetric framework and 
two of the favored supersymmetric WIMP candidates, neutralinos and charged sleptons. 
Following this, we consider WIMPs in extra dimensional scenarios. 



4 



A. Neutralino WIMPs 



A general neutralino x j s a mixture of the neutral Bino, Wino, and Higgsinos. Writing 
x = Nu(-iB) + N 12 (-iW) + N 13 H U + N 1A H d , we find the decay width 

i / v i i i ~ ( u w~ n i i ■ — i- i / v -i i ~ ^ i i i n , i / / / / / ' ~ 
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(8) 



This decay width, and all those that follow, includes the contributions from couplings to 
both the spin ±3/2 and ±1/2 gravitino polarizations. These must all be included, as they 
are comparable in models with high-scale supersymmetry breaking. 

There are also other decay modes. The two-body final states ZG and hG may be kine- 
matically allowed, and three-body final states include HG and qqG. For the WIMP life- 
times we are considering, constraints on electromagnetic energy release from BBN are well- 
studied [13, IE 3 ) but constraints on hadronic cascades are much less certain (22] . Below, 
we assume that electromagnetic cascades are the dominant constraint and provide a careful 
analysis of these bounds. If the hadronic constraint is strong enough to effectively exclude 
two-body decays leading to hadronic energy, our results below are strictly valid only for the 
case x — T) where x — > jG is the only possible two-body decay. If the hadronic constraint is 
strong enough to exclude even three-body hadronic decays, such as 7 — > qqG, the entire neu- 
tralino superWIMP scenario may be excluded, leaving only slepton superWIMP scenarios 
(discussed below) as a viable possibility. Detailed studies of BBN constraints on hadronic 
cascades at r ~ 10 6 s may therefore have important implications for super WIMPs. 

With the above caveats in mind, we now focus on Bino-like neutralinos, the lightest 
neutralinos in many simple supergravity models. For pure Binos, 
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(9) 



In the limit Am = m WIMP — mswiMP <^ ^swimp, — > jG) oc (Am) 3 and the decay 
lifetime is 

100 GeVl 3 



-(B -> 7G) w 2.3 x 10 7 s 



Am 



(10) 



independent of the overall mwiMP, ^swimp mass scale. This threshold behavior, sometimes 
misleadingly described as P-wave, follows not from angular momentum conservation, but 
rather from the fact that the gravitino coupling is dimensional. For the case S = 7, clearly 
all of the initial photon energy is deposited in an electromagnetic shower, so 

e EM = £ 7 , £ had ~ . (11) 

If the WIMP is a Bino, given values of m WIM p and m SWIMP , r is determined by Eq. Q, 
and Eqs. © and determine the energy release Cem- These physical quantities are given 
in Fig. H] for a range of (ms W i M p, Am). 



B. Charged Slepton WIMPs 



For a slepton NLSP, the decay width is 




10 6 10 9 10 12 10 15 10 18 10 6 10 9 10 12 10 15 10 18 
t wimp (sec) t wimp (sec) 



FIG. 1: Predicted values of WIMP lifetime r and electromagnetic energy release £em = £em^wimp 
in the B (left) and f (right) WIMP scenarios for m SW iMP = 1 GeV, 10 GeV, . . . , 100 TeV (top to 
bottom) and Am = mwiMP — w-swimp = 1 TeV, 100 GeV, . . . , 100 MeV (left to right). For the f 
WIMP scenario, we assume £em = \E t . 



This expression is valid for any scalar superpartner decaying to a nearly massless standard 
model partner. In particular, it holds for i — e, ft, or f, and arbitrary mixtures of the £l and 
£r gauge eigenstates. In the limit Am = mwiMP — ^swimp <C ^swimp, the decay lifetime is 



3.6 x 10 s 



100 GeV 



Am 



m G 

1 TeV 



(13) 



For selectrons, the daughter electron will immediately initiate an electromagnetic cascade, 



so 



e E M ^ E e , e had ~ . 



(14) 



Smuons produce muons. For the muon energies E„ ~ M wea k and temperatures T T of interest, 



E^T T -C m . These muons therefore interact with background photons through /x7bg 
with the Thomson cross section for muons. The interaction time is 



Tint — 



avn: 



BG 

7 



'8vra 2 \ /2C(3)T r 3S 



3m2 



7T Z 



~ 7 x 10" 5 



keV 

~t7 



HI 



(15) 



This is typically shorter than the time-dilated muon decay time (E^/m^) 2.0 x 10~ 6 s. The 
muon energy is, therefore, primarily transferred to electromagnetic cascades, and so 



£EM — E/j, , £ n ad — 



(16) 



If muons decay before interacting, some electromagnetic energy will be lost to neutrinos, 
but in any case, e had 0, and hadronic cascades may be safely ignored. 
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Finally, stau NLSPs decay to taus. Before interacting, these decay to e, /x, 7r°, ti ± and v 
decay products. All of the energy carried by e, /i, and n° becomes electromagnetic energy. 
Decays n + — > ji + u also initiate electromagnetic cascades with energy ~ E n +/2. Making the 
crude assumption that energy is divided equally among the r decay products in each decay 
mode, and summing the e, /i, tt°, and half of the 7r ± energies weighted by the appropriate 
branching ratios, we find that the minimum electromagnetic energy produced in r decays 
is ~ \E T . The actual electromagnetic energy may be larger. For example, for charged 
pions, following the analysis for muons above, the interaction time for 7t ± 7bg - ► 7i" ± 7 is of the 
same order as the time-dilated decay time (E^±/m w ±) 2.6 x 10~ 8 s. Which process dominates 
depends on model parameters. Neutrinos may also initiate electromagnetic showers if the 
rate for i^bg ~~ > e + e~ is significant relative to vv-bg vv - 

All of the t decay products decay or interact electromagnetically before initiating hadronic 
cascades. The hadronic interaction time for pions and kaons is 



a"hadf?7n^ G 



18 s 



100 mb 



6 x 10 
V 



-10 



keV 



3 



(17) 
(18) 



where r\ is the baryon-to-photon ratio, and we have normalized the cross section to the 
largest possible value. We see that hadronic interactions are completely negligible, as there 
are very few nucleons with which to interact. In fact, the leading contribution to hadronic 
activity comes not from interactions with the existing baryons, but from decays to three- 
body and four-body final states, such as iZG and iqqG, that may contribute to hadronic 
energy. However, the branching ratios for such decays are also extremely suppressed, with 
values ~ 10~ 3 — 10~ 5 23||. In contrast to the case for neutralinos, then, the constraints on 



electromagnetic energy release are guaranteed to be the most stringent, and constraints on 
hadronic energy release may be safely ignored for slepton WIMP scenarios. 
Combining all of these results for stau NLSPs, we find that 

£em ~ ^ Et ~ Et ' £ had = 0, (19) 

where the range in eem results from the possible variation in electromagnetic energy from 
^ and v decay products. The precise value of £em is in principle calculable once the 
stau's chirality and mass, and the superWIMP mass, are specified. However, as the possible 
variation in £em is not great relative to other effects, we will simply present results below 
for the representative value of £em = \E T . 

For slepton WIMP scenarios, Eq. (jl2j) determines the WIMP lifetime r in terms of m^iMP 
and mswiMP, and Cem is determined by Eq. © and either Eq. (JHJ), (JUJ), or (fTTJj). These 
physical quantities are given in Fig. [T] in the f WIMP scenario for a range of (mwiMP, Am). 
For natural weak-scale values of these parameters, the lifetimes and energy releases in the 
neutralino and stau scenarios are similar. A significant difference is that larger WIMP 
masses are typically required in the slepton scenario to achieve the required relic abundance. 
However, thermal relic densities rely on additional supersymmetry parameters, and such 
model-dependent analyses are beyond the scope of this work. 
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C. KK gauge boson and KK charged lepton WIMPs 

In scenarios with TeV _1 -size universal extra dimensions, KK gravitons are superWIMP 
candidates. The WIMPs that decay to graviton superWIMPs then include the 1st level KK 
partners of gauge bosons and leptons. 

For the KK gauge boson WIMP scenario, letting V 1 = B 1 cos 6^ + W 1 sin 6^, 

cos 2 9 W cos 2 9^ + sin 2 9 W sin 2 Q]^ 



TiV 1 -> 7G 1 ) 



72nM? 



x ■ 



m 



G 1 



1 - 



m 



m 



v 1 



G 1 



m 



v 1 



m 



For a Belike WIMP, this reduces to 



r( £i^ 7G i)_ cos2 ^ m ^ 



1 - 



m 



G 1 



1 + 3 



m 



G 1 



m 2 Bl 



+ 6 



72ttM% m%i [ m\ 
In the limit Am = m WIMP — m SW iMP *C m SWIM p, the decay lifetime is 

100 GeV l3 



v 1 



m 4 Bl 



t(B 1 ^7G 1 ) « 1.4 x 10 7 s 



Am 



(20) 



(21) 



(22) 



independent of the overall mwiMP, "^swimp mass scale, as in the analogous supersymmetric 
case. 

For KK leptons, we have 

4 r 



r^ 1 -> eg 1 ) 



1 



48?rM 2 m^i 



1 



m 



2 + 3 



m 



G 1 



m 2 ! 



(23) 



valid for any KK lepton (or any KK fermion decaying to a massless standard model particle, 
for that matter). In the limit Am = m WIMP — mswiMP <C ^swimp, the decay lifetime is 

T (^ (G ')»r.3xio' s rlooGeV1 ' 



Am 



m G i 
1 TeV 



(24) 



In all cases, the expressions for e EM and £had are identical to those in the analogous super- 
symmetric scenario. 

KK graviton superWIMPs are therefore qualitatively similar to gravitino superWIMPs. 
The expressions for WIMP lifetimes and abundances are similar, differing numerically only 
by 0(1) factors. We therefore concentrate on the supersymmetric scenarios in the rest of 
this paper, with the understanding that all results apply, with 0(1) adjustments, to the case 
of universal extra dimensions. A more important difference is that the desired thermal relic 
density is generally achieved for higher mass WIMPs in extra dimensional scenarios that in 
the supersymmetric case. 

III. BARYOMETRY 

A. Standard BBN and CMB Baryometry 



Big Bang nucleosynthesis predicts primordial light element abundances in terms of one 
free parameter, the baryon-to-photon ratio 77 = n B /n 7 . At present, the observed D, 4 He, 
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3 He, and 7 Li abundances may be accommodated for baryon-to-photon ratios in the range j2- 



7] 10 = 77/IO 



10 



2.6 - 6.2 . 



(25) 



In light of the difficulty of making precise theoretical predictions and reducing (or even 
estimating) systematic uncertainties in the observations, this consistency is a well-known 
triumph of standard Big Bang cosmology. 

At the same time, given recent and expected advances in precision cosmology, the stan- 
dard BBN picture merits close scrutiny. Recently, BBN baryometry has been supplemented 
by CMB data, which alone yields r] w = 6.1 ± 0.4 [9]. Observations of deuterium ab- 
sorption features in spectra from high redshift quasars imply a primordial D fraction of 
D/H = 2.78t° f 8 x 10- 5 @. Combined with standard BBN calculations |U, this yields 
7710 = 5.9 ± 0.5. The remarkable agreement between CMB and D baryometers has two 
new implications for scenarios with late-decaying particles. First, assuming there is no fine- 
tuned cancellation of unrelated effects, it prohibits significant entropy production between 
the times of BBN and decoupling. In Sec. I1V1 we will show that the entropy produced in 
superWIMP decays is indeed negligible. Second, the CMB measurement supports determi- 
nations of 77 from D, already considered by many to be the most reliable BBN baryometer. It 
suggests that if D and another BBN baryometer disagree, the "problem" lies with the other 
light element abundance — either its systematic uncertainties have been underestimated, 
or its value is modified by new astrophysics or particle physics. Such disagreements may 
therefore provide specific evidence for late-decaying particles in general, and superWIMP 
dark matter in particular. We address this possibility here. 

In standard BBN, the baryon-to-photon ratio 7710 = 6.0 ± 0.5 favored by D and CMB 
observations predicts [26| 



prediction of Eq. (J28|) . The 7 Li fraction may be determined precisely in very low metallicity 
stars. Three independent studies find 



where depletion effects have been estimated and included in the last value. Within the 
published uncertainties, the observations are consistent with each other but inconsistent 
with Eq. (|28|) . with central values lower than predicted by a factor of 3 — 4. 7 Li may be 
depleted from its primordial value by astrophysical effects, for example, by rotational mixing 
in stars that brings Lithium to the core where it may be burned [3^,0], but it is controversial 
whether this effect is large enough to reconcile observations with the BBN prediction |29| . 

The other light element abundances are in better agreement. For example, for 4 He, Olive, 
Skillman, and Steigman find Y p = 0.234±0.002 j32|, lower than Eq. (|26p. but the uncertainty 
here is only statistical. Y p is relatively insensitive to 77 and a subsequent study of Izotov and 
Thuan finds the significantly higher range 0.244 ±0.002 (33[. 3 He has recently been restricted 





(29) 
(30) 
(31) 
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to the range 3 He/H < (1.1 ± 0.2) x 10" 5 Q, consistent with the CMB + D prediction of 
Eq. (|27|1 . Given these considerations, we view disagreements in 4 He and 3 He to be absent or 
less worrisome than in 7 Li. This view is supported by the global analysis of Ref. [261 ] . which, 
taking the "high" Y p values of Izotov and Thuan, finds x 2 = 23.2 for 3 degrees of freedom, 
where \ 2 is completely dominated by the 7 Li discrepancy. 



B. Super WIMPs and the 'Li Underabundance 

Given the overall success of BBN, the first implication for new physics is that it should 
not drastically alter any of the light element abundances. This requirement restricts the 
amount of energy released at various times in the history of the universe. A recent analysis 
by Cyburt, Ellis, Fields, and Olive of electromagnetic cascades finds that the shaded regions 
of Fig. El are excluded by such considerations jl7(. The various regions are disfavored by the 
following conservative criteria: 



D low 
D high 
4 He low 
7 Li low 



D/H < 1.3 x 10~ 5 (32) 

D/H > 5.3 x 10~ 5 (33) 

Y p < 0.227 (34) 

7 Li/H < 0.9 x 10~ 10 . (35) 



A subset of superWIMP predictions from Fig. ^ is superimposed on this plot. The 
subset is for weak-scale mgwiMP and Am, the most natural values, given the independent 
motivations for new physics at the weak scale. As discussed previously [If, the BBN con- 
straint eliminates some of the region predicted by the superWIMP scenario, but regions with 
^wimp,^swimp ~ M weak remain viable. 

The 7 Li anomaly discussed above may be taken as evidence for new physics, however. 
To improve the agreement of observations and BBN predictions, it is necessary to destroy 
7 Li without harming the concordance between CMB and other BBN determinations of r\. 
This may be accomplished for (r, (em) ~ (3 x 10 6 s, 10~ 9 GeV), as noted in Ref. [TtJ. This 
"best fit" point is marked in Fig. [21 The amount of energy release is determined by the 
requirement that 7 Li be reduced to observed levels without being completely destroyed - 
one cannot therefore be too far from the " 7 Li low" region. In addition, one cannot destroy 
or create too much of the other elements. 4 He, with a binding threshold energy of 19.8 
MeV, much higher than Lithium's 2.5 MeV, is not significantly destroyed. On the other 
hand, D is loosely bound, with a binding energy of 2.2 MeV. The two primary reactions are 
D destruction through 7D — + np and D creation through 7 4 He — > DD. These are balanced 
in the channel of Fig. [21 between the "low D" and "high D" regions, and the requirement 
that the electromagnetic energy that destroys 7 Li not disturb the D abundance specifies the 
preferred decay time r ~ 3 x 10 6 s. 

Without theoretical guidance, this scenario for resolving the 7 Li abundance is rather 
fine-tuned: possible decay times and energy releases span tens of orders of magnitude, 
and there is no motivation for the specific range of parameters required to resolve BBN 
discrepancies. In the superWIMP scenario, however, both r and (em are specified: the 
decay time is necessarily that of a gravitational decay of a weak-scale mass particle, leading 
to Eq. (Q), and the energy release is determined by the requirement that super WIMPs be 
the dark matter, leading to Eq. (J7|). Remarkably, these values coincide with the best fit 
values for r and (em- More quantitatively, we note that the grids of predictions for the B 
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FIG. 2: The grid gives predicted values of WIMP lifetime r and electromagnetic energy release 
(em = eem^wimp in the B (left) and f (right) WIMP scenarios for mswiMP = 100 GeV, 300 GeV, 
500 GeV, 1 TeV, and 3 TeV (top to bottom) and Am = m wm p - ^swimp = 600 GeV, 400 GeV, 
200 GeV, and 100 GeV (left to right). For the f WIMP scenario, we assume e E M = \E T . The 
analysis of BBN constraints by Cyburt, Ellis, Fields, and Olive jl^J excludes the shaded regions. 
The best fit region with (r, Cem) ~ (3 x 10 6 s, 10 -9 GeV), where 7 Li is reduced to observed levels 
by late decays of WIMPs to superWIMPs, is given by the circle. 



and f scenarios given in Fig. |21 cover the best fit region. Current discrepancies in BBN light 
element abundances may therefore be naturally explained by superWIMP dark matter. 

This tentative evidence may be reinforced or disfavored in a number of ways. Improve- 
ments in the BBN observations discussed above may show if the 7 Li abundance is truly below 
predictions. In addition, measurements of 6 Li/H and 6 Li/ 7 Li may constrain astrophysical 
depletion of 7 Li and may also provide additional evidence for late decaying particles in 
the best fit region 

Finally, if the best fit region is indeed realized by 



WIMP — > SWIMP decays, there are a number of other testable implications for cosmology 
and particle physics. We discuss these in Sees. |V] and EU 



IV. ENTROPY PRODUCTION 



In principle, there is no reason for the BBN and CMB determinations of r] to agree - 
they measure the same quantity, but at different epochs in the universe's history, and rj may 
vary [21]. What is expected to be constant is the number of baryons 



N 



B 



n B R 3 



2C(3) 



T i Ri 



7T Z 



(36) 
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where R is the scale factor of the universe. Since the entropy S is proportional to T 3 R 3 
when g ts , the number of relativistic degrees of freedom for entropy, is constant, 

^ = f , (37) 

where the superscripts and subscripts i and / denote quantities at times ti and tf, respec- 
tively. The quantities rji and r\f therefore must agree only if there is no entropy production 
between times U and tf. 

Conversely, as noted in Sec. the agreement of CMB and D baryometers implies 
that there cannot be large entropy generation in the intervening times |2l| . barring fine- 
tuned cancellations between this and other effects. WIMP decays occur between BBN and 
decoupling and produce entropy. In this section, we show that, for energy releases allowed 
by the BBN constraints discussed above, the entropy generation has a negligible effect on 
baryometry. 

We would like to determine the change in entropy from BBN at time ti to decoupling at 
time tf. The differential change in entropy in a comoving volume at temperature T is 

dS = ^ , (38) 
where the differential energy injected into radiation is 

dQ = £em^wimp-R 3 — • (39) 

T 

In Eq. (I39|) . n-wiMP is the WIMP number density per comoving volume. R may be eliminated 
using 

2vr 2 

S = -^9*sT 3 R 3 . (40) 
Substituting Eqs. and (j4Ti|) into Eq. (J3HJ) and integrating, we find 

(41) 



Sf 

— = exp 



■f 45 1 dt 

^EM^WIMP ■ 



ti 



27T T A T 



As WIMPs decay, their number density is 



p3 „ q 7^3 

"-WIMP = n WIMP^ e = n WIMP ~i ^3 e ) l 4 ^J 



and so 

s~ = exp 



45 1 ftt S t T { _ t/T dt 
£EMnwiMP 2^gl s TfJu ST 6 r 



(43) 



Equation (}4*3"|) is always valid. However, it is particularly useful if the change in entropy 
may be treated as a perturbation, with AS S^. Given the high level of consistency of 77 
measurements from deuterium and the CMB, this is now a perfectly reasonable assumption. 
We may therefore solve Eq. (|43|) iteratively. In fact, the first approximate solution, obtained 
by setting Si/ S = 1 in the integral, is already quite accurate. The integral may be further 
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simplified if the universe is always radiation dominated between BBN and decoupling. This 
is certainly true in the present analysis, as 



Pwimp v 60 C (3) 

- mwiMP i wimp - 



PR 



m WIMP^WIMP 



4.5 x 10~ 6 GeV 



3.36 

9* 



1 keV 
T 



< 1 . 



(44) 



WIMPs therefore decay before their matter density dominates the energy density of the 
universe. We may then use the radiation-dominated era relations 



t 



1 
2H 



H 2 



8tt 



PR 



PR 



30 



(45) 



to eliminate T in favor of t in the integral of Eq. (}4*3*|) . Finally, < r < t/, and, as the 
dominant contribution to the integral is from t ~ r, we may replace g* by gl, its (constant) 
value during the era of WIMP decay. 

Exploiting all of these simplifications, the integral in Eq. (J43j) reduces to 



U9*\' A f±_) l/2 e -t/rdt 



T 




t\^ e _ t/T dt 
T 



(46) 
(47) 



Finally, substituting Eq. (}4Tj) into Eq. 
relations of Eq. (|45jh we find 



Si 

For small entropy changes, 
AS 



cxp 



C(3) 



and again using the radiation-dominated era 

45 3 / 4 a?n i/4 £EMnWiMP 



71 



11/4 



9ls 



Mi 



pi 



S; 



ln s 



1.10 x 10" 



EM 



10" 9 GeV 



10 6 s 



(48) 



(49) 



where we have used C(3) ~ 1.202, and gl ~ 3.36 and g* s ~ 3.91 are the appropriate degrees 
of freedom, which include only the photon and neutrinos. 

Contours of AS/ Si are given in the (r, Cem) plane in Fig. El for late-decaying Binos and 
staus. For reference, the BBN excluded and best fit regions are also repeated from Fig. as 
are the regions predicted for natural superWIMP scenarios. We find that the superWIMP 
scenario naturally predicts AS/ Si ^$ 10~ 3 . Such deviations are beyond foreseeable sensi- 
tivities in studies CMB and BBN baryometry. Within achievable precisions, then, CMB 
and BBN baryometers may be directly compared to each other in superWIMP dark matter 
discussions, as we have already done in Sec. IIHI 

Entropy production at the percent level may be accessible in future baryometry studies. 
It is noteworthy, however, that, independent of theoretical framework, such large entropy 
production from electromagnetic energy release in late-decaying particles is excluded by 
BBN constraints for decay times 10 4 s < r < 10 12 s. Only for decays very soon after BBN 
times ti ~ 1 — 100 s or just before decoupling times tf ~ 10 13 s can entropy production 
significantly distort the comparison between BBN and CMB baryon-to-photon ratios. In 
fact, only the very early decays are a viable source of entropy production, as very late time 
decays create unobserved CMB black body distortions, which we now discuss. 
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FIG. 3: Contours of fractional entropy production AS/ Si from late decays in the (r, Cem) plane. 
Regions predicted by the superWIMP dark matter scenario and BBN excluded and best fit regions 
are given as in Fig. |2J 



V. IMPLICATIONS FOR CMB BLACK BODY DISTORTIONS 

The injection of electromagnetic energy may also distort the frequency dependence of the 
CMB black body radiation. For the decay times of interest, with redshifts z ~ 10 5 — 10 7 , the 
resulting photons interact efficiently through — > 7e _ , but photon number is conserved, 
since double Compton scattering — > 77c - and thermal bremsstrahlung eX — > eX'y, 
where X is an ion, are inefficient. The spectrum therefore relaxes to statistical but not 
thermodynamic equilibrium, resulting in a Bose-Einstein distribution function 



e E/(kT)+tM 



(50) 



with chemical potential /1 7^ 0. 

For the low values of baryon density currently favored, the effects of double Compton 
scattering are more significant than those of thermal bremsstrahlung. The value of the 
chemical potential \x may therefore be approximated for small energy releases by the analytic 
expression 



H = 8.0 x 10- 



where 



In Fif 

10~ 5 E 



r dC = 6.1 x 



10 6 s 



T 


1 

2 


Cem 




-(rac/r) 5 / 4 


LlO 6 s. 




10- 


9 GeV_ 


e 


r T 




_ 12 
5 


~tt B h 2 ~ 


4 
5 


n - iy 1 

2 P 


.2.725 K. 




0.022 




0.88 



(51) 



(52) 



we show contours of chemical potential \x. The current bound is fi < 9 x 
. We see that, although there are at present no indications of deviations from 
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FIG. 4: Contours of jx, parameterizing the distortion of the CMB from a Planckian spectrum, in 
the (r, Cem) plane. Regions predicted by the superWIMP dark matter scenario, and BBN excluded 
and best fit regions are given as in Fig. [2J 

black body, current limits are already sensitive to the superWIMP scenario, and particularly 
to regions favored by the BBN considerations described in Sec. 11111 In the future, the Diffuse 
Microwave Emission Survey (DIMES) may improve sensitivities to /x ~ 2 x 10 -6 |38|. DIMES 
will therefore probe further into superWIMP parameter space, and will effectively probe all 
of the favored region where the 7 Li underabundance is explained by decays to superWIMPs. 

VI. IMPLICATIONS FOR PARTICLE PHYSICS 

The superWIMP scenario has implications for the superpartner (and KK) spectrum, and 
for searches for supersymmetry (and extra dimensions) at particle physics experiments. In 
this section, we consider some of the implications for high energy colliders. 

Lifetimes and energy releases are given as functions of mswiMP and Am in Fig. BBN 
and CMB baryometry, along with limits on CMB [i distortions, exclude some of this param- 
eter space. The excluded regions were presented and discussed in Ref. 

Here we concentrate on the regions preferred by the tentative evidence for late decaying 
particles from BBN considerations. As noted above, the preferred lifetimes and energy 
releases for which 7 Li is reduced without sacrificing the concordance between CMB and D r\ 
determinations is a region around (r, Cem) ~ (3 x 10 6 s, 10~ 9 GeV). This region is highlighted 
in Fig. For the f case, we present a range of best fit regions to account for the possible 
range £em = (§ — 1)E T of Eq. (fTTJj) discussed in Sec. HU 

Given some variation in the preferred values of r and Cem, there is a fair amount of 
variation in the underlying superpartner masses. We may draw some rough conclusions, 
however. For the B WIMP scenario the preferred parameters are ttiq ~ 600 GeV and 
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FIG. 5: Contours of constant r (dashed, red) and constant Cem = ^em^wimp (solid, blue) in 
the (mswiMP > Am) plane in the B (left) and f (right) WIMP scenarios. The regions with BBN 
preferred values (r, Cem) ~ (3 x 10 6 s, 10 -9 GeV) are given by the circles. For the f WIMP scenario, 
the solid circle is favored if £em = \E T \ the dashed circles are favored if £em = \E T or £em = E T . 

trig ~ 800 GeV. All other superpartners are necessarily heavier than m^. The resulting 
superpartner spectrum is fairly heavy, although well within reach of the LHC, assuming the 
remaining superpartners are not much heavier. This scenario will be indistinguishable at 
colliders from the usual supergravity scenario where the gravitino is heavier than the LSP 
and the usual signal of missing energy from neutralinos applies. 

For the f superWIMP scenario, there are dramatic differences. From Fig. El the BBN 
preferred masses are tuq ~ 300 — 1100 GeV and Am = — ~ 300 — 400 GeV. Although 
fairly heavy, this range of superpartner masses is again well within the reach of the LHC 
and possibly even future linear colliders. In this case, collider signatures contrast sharply 
with those of standard supergravity scenarios. Typically, the region of parameter space in 
which a stau is the lightest standard model superpartner is considered excluded by searches 
for charged dark matter. In the superWIMP scenario, this region is allowed, as the stau is 
not stable, but metastable. Such particles therefore evade cosmological constraints, but are 
effectively stable on collider time scales. They appear as slow, highly-ionizing charged tracks. 
This spectacular signal has been studied in the context of gauge-mediated supersymmetry 
breaking models with a relatively high supersymmetry-breaking scale [39], and discovery 
limits are, not surprisingly, much higher than in standard scenarios. It would be interesting 
to evaluate the prospects for discovering and studying meta-stable staus at the Tevatron, 
LHC, and future linear colliders in various superWIMP scenarios. 
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VII. CONCLUSIONS AND FUTURE DIRECTIONS 



SuperWIMP dark matter presents a qualitatively new dark matter possibility realized 
in some of the most promising frameworks for new physics. In supergravity, for example, 
superWIMP dark matter is realized simply by assuming that the gravitino is the LSP. When 
the NLSP is a weakly-interacting superpartner, the gravitino superWIMP naturally inherits 
the desired dark matter relic density. The prime WIMP virtue connecting weak scale physics 
with the observed dark matter density is therefore preserved by superWIMP dark matter. 

Because superWIMP dark matter interacts only gravitationally, searches for its effects 
in standard dark matter experiments are hopeless. At the same time, this superweak in- 
teraction implies that WIMPs decaying to it do so after BBN. BBN observations and later 
observations, such as of the CMB, therefore bracket the era of WIMP decays, and provide 
new signals. SuperWIMP and conventional WIMP dark matter therefore have disjoint sets 
of signatures, and we have explored the new opportunities presented by superWIMPs in this 
study. We find that the superWIMP scenario is not far beyond reach. In fact, precision 
cosmology already excludes some of the natural parameter space, and future improvements 
in BBN baryometry and probes of CMB \i distortions will extend this sensitivity. 

We have also found that the decay times and energy releases generic in the superWIMP 
scenario may naturally reduce 7 Li abundances to the observed levels without sacrificing the 
agreement between D and CMB baryometry. The currently observed 7 Li underabundance 
therefore provides evidence for the superWIMP hypothesis. This scenario predicts that more 
precise BBN observations will expose a truly physical underabundance of 7 Li. In addition, 
probes of CMB fi distortions at the level of fi ~ 2 x 10~ 6 will be sensitive to the entire 
preferred region. An absence of such effects will exclude this explanation. 

We have considered here the cases where neutralinos and sleptons decay to gravitinos and 
electromagnetic energy. In the case of selectrons, smuons, and staus, we have shown that 
BBN constraints on electromagnetic cascades provide the dominant bound. For neutralinos, 
however, the case is less clear. Neutralinos may produce hadronic energy through two-body 
decays x ~^ ZG, hG, and three-body decays x qqG. Detailed BBN studies constraining 
hadronic energy release may exclude such two-body decays, thereby limiting possible neu- 
tralino WIMP candidates to photinos, or even exclude three-body decays, thereby eliminat- 
ing the neutralino WIMP scenario altogether. At present, detailed BBN studies of hadronic 
energy release incorporating the latest data are limited to decay times r < 10 4 s j22|. We 
strongly encourage detailed studies for later times r ~ 10 6 s, as these may have a great 
impact on what superWIMP scenarios are viable. 

Finally, in the course of this study, we presented a model-independent study of entropy 
production in light of the recent WMAP data. The agreement of precise CMB and D 
baryon-to-photon ratios limits entropy production in the time between BBN and decou- 
pling. However, constraints on BBN light element abundances and CMB distortions already 
provide stringent bounds. We have compared these constraints here. We find that BBN 
abundances and CMB black body distortions largely eliminate the possibility of significant 
entropy production. For fractional entropy changes at the percent level, which may be visible 
through comparison of future BBN and CMB baryometers, these other constraints require 
the entropy production to take place before ~ 10 4 s, that is, in a narrow window not long 
after BBN. 
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